Separation of microscale chiral objects by shear flow 
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We show that plane parabolic flow in a microfluidic channel causes nonmotile helically-shaped 
bacteria to drift perpendicular to the shear plane. Net drift results from the preferential alignment 
of helices with streamlines, with a direction that depends on the chirality of the helix and the sign of 
the shear rate. The drift is in good agreement with a model based on resistive force theory, and sep- 
aration is efficient (> 80%) and fast (< 2 s). We estimate the effect of Brownian rotational diffusion 
on chiral separation and show how this method can be extended to separate chiral molecules. 



Many biochemically active molecules are naturally chi- 
ral and can only bind to target chiral molecules of a 
specific handedness The other enantiomer (i.e. the 
molecule having opposite handedness) may be inactive 
or cause undesirable effects. Chemical synthesis of chiral 
molecules usually produces a racemic mixture, with equal 
amounts of both enantiomers, and their separation based 
on chirality is of importance in fields ranging from agri- 
culture to food and pharmaceutical industries. Currently 
favored approaches rely on gas, liquid or capillary elec- 
tromigration chromatography [2], requiring costly chiral 
media. Thus, simpler, alternative approaches to chiral 
separation are desirable. 

Several alternative proposals for chiral separation ex- 
ploit hydrodynamic forces. Some of these, yet untested 
experimentally, rely on the presence of a surface [3| or 
array of microvortices [4|, and there has been success- 
ful chiral separation of cm-sized crystals in a rotating 
drum [5]. Other methods [6, 7] stem from the prediction 
that a chiral particle in a simple shear flow experiences a 
lateral drift [8j- However, the feasibility of this approach 
has remained questionable, as measurements in Couette 
cells reported that the drift of mm-sized chiral objects [si 
and the forces on cm-sized ones [To[ differ from predic- 
tions by two orders of magnitude [q| or even in sign [To| . 

Here we report that microscale chiral objects, three or- 
ders of magnitude smaller than previous studies 0, [Toj| , 
experience a lateral drift in a microfluidic shear flow and 
the magnitude of the drift is in agreement with our the- 
ory. Previous work has demonstrated the ability of mi- 
crofluidics to separate and sort colloids by size , sper- 
matozoa by motility [12] . and microbes by their prefer- 
ence for dissolved chemicals [13|- Our method uses mi- 
crochannels to sort particles by chirality. We show that 
an enantiomer drifts with direction determined by the lo- 
cal shear, demonstrate the feasibility of this method for 
chiral separation, and indicate how the high shear rates 
achievable in microchannels (> 10 6 s _1 |l4|) allow it to 
be extended to smaller scales (< 40 nm). 

The origin of chirality-dependent drift at low Reynolds 
number can be simply understood for the case of a he- 
lix. In a shear flow, objects undergo periodic rotations 




FIG. 1: (color) Schematic of a right-handed helix in simple 
shear flow. Red and black shadings show the top and bottom 
halves of the helix. Upper inset: the net force acting on 
one pitch of the helix is along — z. Lower inset: predicted 
normalized drift velocity v/vo vs. helix orientation 6. 6 is the 
angle between a helix in the x-y plane and the flow, and vo is 
the drift velocity of a helix aligned with the flow (0 = 0). 



known as Jeffery orbits [l5| : a sphere rotates with con- 
stant angular velocity, whereas for an elongated body, 
such as a helix, the velocity depends on orientation. The 
more elongated a body, the longer its residence time when 
aligned with streamlines. In addition to rotating in a Jef- 
fery orbit, a helix drifts across streamlines. To see why, 
consider a right-handed helix aligned with a simple shear 
flow (Fig.[T]), and decompose the velocity at a segment of 
the helix into components perpendicular and parallel to 
the segment. Drag on a thin rod in low Reynolds number 
flow is anisotropic, with a greater resistance when ori- 
ented perpendicular rather than parallel to the flow [l6| . 
Since the flows at the top (red) and bottom (black) halves 
of the helix are in opposite directions, both halves have 
a drag component along — z (Fig. [TJ top inset). Thus, a 
right-handed helix aligned with the shear flow drifts in 
the — z direction. Reversing the chirality of the helix or 
the sign of the shear produces a drift along +z. Further- 
more, the drift depends on the orientation of the helix. 
A right-handed helix aligned with the y axis drifts along 
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FIG. 2: (color) (a) Scanning Electron Micrograph (SEM) of L. biflexa flaB mutant, with typical dimensions (inset). The 
bent configuration is a result of SEM preparation and live organisms are nearly always straight, (b) MicroChannel design with 
separate inlets for spirochetes and buffer. The color-coded squares refer to the locations of data collection (Fig. 3b). (c) 
Schematic of the separation process (for one enantiomer) in the microchannel (W — 1mm, H — 90 /im, w « 100 /im). The 
lateral drift direction depends on the sign of the shear 7, resulting in divergence of top and bottom streams. 



+z. In general, the drift velocity v of a helix lying in the 
x-y plane is a function of the angle between the helix 
and the flow (Fig. [TJ bottom inset). Hence, the mean 
drift v depends on the distribution of the orientation of 
the helix relative to the flow. When the helix is pref- 
erentially aligned with streamlines, for example due to 
JefTery orbits, it has a net drift. When all orientations 
are equally likely, for example due to strong Brownian 
rotational diffusion, the net drift is zero [7]. Here we use 
microscale objects to demonstrate the chirality-induced 
drift and determine the size of the smallest helices that 
can be separated using this principle. 

We use a nonmotile, right-handed, helically-shaped 
bacterium, Leptospira biflexa flaB mutant [17] (Fig. 2a), 
as the microscale chiral object. These spirochetes have 
an average length of 16 /im, thickness of 150 nm, and di- 
ameter of 200 nm (Fig. [2^ and [18[). To expose spiro- 
chetes to a shear flow, we fabricated a 110 cm long, 
serpentine-shaped microfluidic channel via soft lithogra- 
phy [l9|| (Figs. [2)3, c). The channel has rectangular cross 
section of depth H = 90 /im and width W = 1 mm, re- 
sulting in a parabolic flow, u(y) = (3/7/2) [1 — (2y/H) 2 ], 
with mean speed U = 3.09 mms -1 . Due to the high as- 
pect ratio W/H, the velocity u is uniform in z, except 
for a 150 /im layer adjacent to the two sidewalls [20|, 
inconsequential for these experiments. The shear rate 
j(y) = du/dy = —12yU/H 2 varies linearly along y (be- 
tween ±206 s _1 ), with 7 < for y > 0, and 7 > for 
y < (Fig. [2fc). The microchannel is equipped with 
two inlets, one for spirochetes and one for buffer solu- 
tion (EMJH liquid medium [I?})- Spirochetes were intro- 
duced through a microinjector, forming a band of width 
w w 100 /im at the center of the channel, uniform over 
the depth, with buffer streams on either side. As 7 varies 
with one expects the drift to vary with vertical posi- 
tion, and spirochetes in the top and bottom halves of the 
channel to drift in opposite directions. 

We imaged the spirochete population using phase con- 



trast microscopy with a 40 x objective and a CCD cam- 
era. At ten locations along the channel (Fig. Od), sets 
of fifty images were recorded at 1.5 s intervals at three 
depths: y = -H/4 ("bottom"), ("mid"), and H/4 
("top"). Experiments revealed a clear and reproducible 
drift, with direction determined by the sign of the local 
shear rate. A superposition of fifty images acquired at 
x = 107.5 cm shows drift of comparable magnitude and 
opposite direction in the top and bottom halves, and no 
drift at mid-plane, where 7 = (Fig. [3b and supple- 
mental movie [2lj). Using image analysis to obtain the 
z position of individual spirochetes, we determined the 
probability density function (pdf) of the across-channel 
distribution of cells, their mean position z, and the one- 
standard-deviation spread a, for each imaging location 
(x,y). The drift of the top and bottom populations in- 
creases monotonically with x or, equivalent ly, with time. 
The apparent spread of the spirochete distribution with 
increasing x (Fig. [3b, c) is a result of the variation of drift 
with depth combined with the objective's finite depth of 
focus 21]. Brownian motion was negligible (a w 6 /im 
over x = 107.5 cm (2lj). 

Since drift (z) increases linearly with distance x along 
the channel (Fig. 3d), the mean drift velocity is constant 
(v = 0.41 /im s _1 ), resulting in a linearly increasing sep- 
aration distance Az between the top and bottom pop- 
ulations (Fig. [3^). Control experiments with neutrally 
buoyant 1 /im spheres revealed no drift, irrespective of the 
vertical position in the channel (Fig.(3]i, empty squares), 
and thus no separation (Fig.[3k, empty squares), support- 
ing the conclusion that drift is a result of chirality. We 
further verified that the 180-degree turns in the channel 
are inconsequential for these observations (2lj . 

To quantify how the bottom and top streams diverge, 
we define a separation efficiency r\ = \Nl — Nr\/(Nl + 
Nr), where Nl and Nr are the numbers of spirochetes 
to the left and right of the mean position z of the mid 
population (y = 0), and 77 ranges from zero for no separa- 
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FIG. 3: (color) The observed spirochete distribution across the channel at (a) x = 2 cm and (c) x = 107.5 cm. Colors represent 
populations at different depths y: top (orange), mid (black) and bottom (blue) (see Fig. 2c). The mid point of each rectangle 
corresponds to the mean position (z) of the distribution, while the half- width is two standard deviations (2a). (b) Pdf of 
spirochete distribution at various distances x along the channel (see Fig. [2j)). The thick pink and green profiles correspond 
to panels (a) and (c), respectively. (d,e,f) Quantification of lateral drift and separation along x: (d) Mean position of the 
population at three depths, (e) Divergence between the top and bottom populations, (f) Separation efficiency. In all panels, 
full circles refer to spirochete data, lines are linear fits, and empty squares are experiments with spherical beads. Over 10,000 
spirochetes were imaged and located at each position (x,y). 



tion to one for perfect separation. We found r] to increase 
with x for spirochetes, reaching a plateau at r] w 0.8 af- 
ter x = 70 cm, while r] < 0.2 for spheres (Fig. [3J, empty 
squares). This separation was achieved in 227 s, and ex- 
periments with U = 600 mm s _1 (7 between ±4 x 10 4 
s _1 ) further accelerated separation (< 2 s). While we 
have separated a single enantiomer based on the sign of 
the shear, by symmetry our results also apply to separa- 
tion of a racemic mixture in shear of a given sign. 

To confirm that the observed lateral drift is due to 
the coupling of shear and chirality, we use resistive force 
theory (RFT) [l6| and model a spirochete as a 25-turn 
helix with dimensions shown in Fig. [2^. The drag force 
per unit length is f = C|| v || + C± v l? where £y, (± and 
vjj, are the resistive coefficients and relative veloci- 
ties parallel and perpendicular to a given segment of the 
helix, respectively. The drag anisotropy Cl/C|| depends 
on the thickness of the spirochete body and varies be- 
tween 1.4 and 1.7 [l6[. In Stokes flow, the linear and 
angular velocities of the helix are determined by total 
force and total moment balance. This yields the drift 
velocity v of a spirochete for an arbitrary orientation rel- 
ative to the flow. Since we observe that the spirochetes 
lie in the x-y plane to within 10 degrees, as expected 
on theoretical grounds [2l|, we adopt a two-dimensional 
model with orientations limited to the x-y plane. The 
mean drift velocity then depends on the probability dis- 
tribution c(0) of orientations of the helix relative to the 
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FIG. 4: (color online) Size limit for separation of isometric he- 
lices of length L and equivalent aspect ratio r = 70. Dashed 
contours show constant Re (spaced by a factor of 10) and 
the thick contour is Re = 0.1. Solid contours show constant 
Pe. The grayscale shows v/vq, where vo is the drift velocity 
of a helix aligned with the flow. The smallest helices which 
can be separated for a given shear rate, determined by v/vo 
= 0.66 (Pe = 50, marked by circle on inset) and Re < 0.1, 
have L ^ 20 (square), 80 (circle), and 400 nm (diamond) for 
7 w 10 8 , 10 6 , andl0 4 s _1 , respectively. Full symbols mark 
the parameter regimes of our experiments: Re = 0.03 (trian- 
gle) and Re = 5 (hexagon). Inset: v/vo vs. Pe. 



flow. This distribution obeys dtc = DrOqC — do(6c), 
where Dr = 2.3 x 10 _3 s _1 is the rotational diffusiv- 
ity of an L =16 /im long spirochete [21( and is the 
rotation rate caused by the shear flow. We approxi- 



mate 6 as the rotation rate of an ellipsoid with aspect 
ratio r, 6 = -j(cos 2 + r 2 sin 2 (9)/(l + r 2 ) [15], with 
r = 70, as computed from the ratio of the rotation rates 
of the spirochetes at = and = tt/2 obtained from 
RFT. The steady-state solution (dtc = 0) for c(0) was 
obtained numerically and the mean drift velocity com- 
puted as v = J v(0)c(0) dO. This resulted in v = 0.81- 
1.22 /im s _1 , which is of the same order of the measured 
value (0.41 fim s _1 ). The residual discrepancy might be 
associated with the flexibility of the spirochetes and ir- 
regularities in their geometry. 

For chirality-induced separation to take place, prefer- 
ential alignment of helices with the flow is required. In 
our experiments, spirochetes were aligned hydrodynami- 
cally via Jeffery orbits. For smaller particles, rotational 
Brownian motion is more effective at randomizing orien- 
tation. To determine the size of the smallest helical par- 
ticles that can be separated in this manner, we find the 
steady distribution for helices with r = 70 and length 
L. The mean drift velocity v is uniquely determined 
by the Peclet number Pe = j/Dr (Fig. [H inset). For 
Pe ^> 1, alignment by shear dwarfs rotational diffusion 
and u«^0) where is the drift velocity of a helix aligned 
with the flow. Our experiments are in this limit, with 
Pe w 4 x 10 4 (7 « 103s- 1 at y = H/4). In contrast, 
for Pe < 10, diffusion destroys the alignment and v — >> 0. 
Assuming isometrically scaled helices of length L, Dr 
scales as L~ 3 , decreasing Pe for smaller helices. This can 
be partially counteracted by increasing 7, but the parti- 
cle Reynolds number Re = 7L 2 jv [ll[ must remain small 
(Re <C 1) to operate in the Stokes regime. Hence, as L 
decreases, 7 can be increased as L~ 2 . The competition 
between Pe and Re determines the smallest helices that 
can be separated for a given shear rate, which we esti- 
mate as L « 20, 80, and 400 nm for 7 w 10 8 , 10 6 , and 10 4 
s _1 , respectively (Fig. |4j open symbols). Shear rates as 
high as 10 6 s _1 have been reported [14] in micro channels, 
and current technology allows pressure heads of 2 x 10 8 
Pa (Microfluidics Corp., Newton, MA), such that a 10 
/im deep, 5 cm long channel could generate shear rates 
of 10 T s _1 , sufficient to separate 40 nm particles. 

While we have demonstrated the lateral drift of chi- 
ral objects, one is ultimately interested in separating a 
mixture of enantiomers. Because the top and bottom 
halves of the channel produce drift in opposite directions, 
a racemic mixture would separate into four quadrants, 
with opposite quadrants containing the same enantiomer. 
We are currently implementing three-dimensional micro- 
fabrication to separately collect particles from the four 
quadrants. While we believe our technique represents a 
significant step forward in enabling reliable chiral sep- 
aration by shear, it remains to be seen whether it can 
extend to the molecular scale (< 10 nm). To this end, in 
addition to increasing shear rates via improved manufac- 
turing techniques, it will be interesting to explore the fea- 
sibility of using alternative alignment mechanisms, such 



as external fields or boundary effects [9[, to counteract 
rotational diffusion. Alternatively, one could operate at 
higher Reynolds number: preliminary experiments show 
successful separation at Re = 5, but further research is 
required to elucidate the role of inertia on separation and 
any constraints on microchannel operation, such as Dean 
vortices and performance under high pressure. Finally, 
not all chiral particles are helical, so it would be useful to 
understand what geometries lead to the best separation. 
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